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A complete next-to-leading order QCD description of resonant Z' production 

and decay into tt final states 
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We discuss QCD radiative corrections to the production of a heavy neutral resonance Z' at the 
LHC assuming that it decays into a ti final state. Compared to previous studies, our computation 
includes top quark decays as well as interference between the Z' signal process and the QCD ti 
background. The interference contribution appears for the first time at next-to-leading order QCD 
and requires new one-loop amplitudes that are not present when signal and background are treated 
separately. We describe some examples of how QCD radiative corrections may influence both the 
exclusion limits and studies of properties of the new resonance, once it is discovered. 



I. INTRODUCTION 

Searches for physics beyond the Standard Model 
(BSM) at the LHC are starting to reach into the TeV 
mass range and exclude large fractions of parameter 
spaces of most popular models of new physics. Many 
of these exclusion limits that often involve rather com- 
plicated final states and extreme energies-are based on 
leading order calculations which are not necessarily ad- 
equate in such situations. While, in general, it is unre- 
alistic to expect that theoretical description of collider 
processes that are relevant for all interesting BSM sce- 
narios will be advanced to next-to-leading orded (NLO) 
QCD accuracy, this can be done for such extensions of 
the SM that are considered to be sufficiently robust. In 
fact, such calculations have been reported recently for a 
variety of new physics models in Refs. [IHEj]. 

One such robust scenario is an extension of the Stan- 
dard Model by an additional £7(1) gauge group which 
generically leads to the appearance of a massive neutral 
flavor-conserving spin-one particle which we will refer to 
as Z' . Physics of such particles has been studied in great 
detail, and it is beyond the scope of this paper to dis- 
cuss all possible model-building aspects. The interested 
reader can find studies in this direction in [(l-[lo| and ref- 
erences therein. For recent discussions in the context of 
the LHC see e.g. Refs. 

This new Z' particle will, generically, couple to all 
quarks, including the top quark. For this reason, it can 
be searched for in the process pp —> Z' — > ti at the LHC. 
In fact, searches for Z' resonances in this process are well 
under way [l7l - [22l ]. Currently, such searches exclude Z' 
with masses smaller than 1.5 — 2 TeV where the exact 
value of the exclusion bound depends on the assumed 
couplings of the Z' to light quarks and the top quark. 
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In this paper we study the QCD corrections to the pro- 
duction of a Z' resonance with its subsequent decay to a 
ti pair, without the simplifying assumptions of previous 
studies. In particular, we include NLO QCD corrections 
to decays of top quarks into observable final states, ac- 
count for all spin correlations, and consider an interfer- 
ence of the signal qq — > Z' — > tt with the gluon-mediated 
QCD background qq — > g* — > ti which appears at next- 
to- leading order for the first time. 

Our computation can be used to perform realistic stud- 
ies of both the discovery potential of Z' — > ti at the 
LHC and precision studies of the properties of the Z' 
resonance once it is discovered. Hence our computation 
complements similar precision studies of Z' collider phe- 
nomenology in Drell-Yan final states [lU The po- 
tential applications of our calculation are illustrated by 
computing if-factors for experimental cuts that are sim- 
ilar to the ones employed by ATLAS and CMS \rM2^ in 
their searches for Z' — > ti and by studying the kinematic 
distribution of a suitably defined relative azimuthal an- 
gle of leptons from semileptonic decays of top quarks [23[ 
that can be used to discriminate between vector and axial 
couplings of Z' to top quarks. 

The rest of the paper is organized as follows. In Sec.lTTl 
we describe the Z' model that we use throughout the 
paper and comment on the general setup of the compu- 
tation. In Sec. Mil we discuss aspects of the NLO QCD 
calculation that are relevant for our analysis. In Sec. IIVI 
we present the numerical results. We conclude in Sec. [V] 



II. SETUP 

In this section we describe the setup of the calculation. 
Our calculation can accommodate any possible renormal- 
izable interaction of the Z' to fermions. Furthermore, we 
take both the Z' mass and its width as input parameters 
of the theory, which allows us to describe with NLO accu- 
racy a generic Z' model. To present our results however, 
we restrict ourselves to the leptophobic top-color model 
of Refs. '11, 24], which is the reference model used by the 
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Narrow Z', r^, 0) 


= 15 GeV 


Broad Z' , V 


LO) _ 

Z' 


150 GeV 




(1,1) 


(-1,-1) 


(1,0) 


(0,-1) 


(1,1) 


(-1,-1) 


(1,0) 


(o,-i) 


LO, fj, = m z > /2 


3.94 


3.76 


5.73 


2.13 


39.0 


36.4 


56.6 


20.9 


LO, /i = m Z i 


3.70 


3.53 


5.39 


2.00 


36.8 


34.3 


53.4 


19.6 


LO, /i = 2m Z / 


3.49 


3.32 


5.08 


1.88 


34.7 


32.4 


50.5 


18.5 


NLO, n = m Z '/2 


4.20 


3.90 


6.11 


2.21 


42.0 


38.2 


61.0 


21.9 


NLO, n = m Z i 


4.17 


3.89 


6.06 


2.20 


41.7 


38.1 


60.5 


21.8 


NLO, /n = 2m z , 


4.12 


3.86 


5.98 


2.18 


41.3 


37.8 


59.9 


21.7 


if— factor, fj, = 


1.07 


1.04 


1.07 


1.04 


1.08 


1.05 


1.08 


1.05 


if— factor, /j, — mz' 


1.13 


1.10 


1.12 


1.10 


1.13 


1.11 


1.13 


1.11 


if —factor, = 2mz' 


1.18 


1.16 


1.17 


1.16 


1.19 


1.17 


1.18 


1.17 



TABLE I: Cross sections (in fb) for pp — >• Z' — > tt — >• llvvjj + X at leading and next-to-leading order in perturbative QCD for 
various Z' models at the 14 TeV LHC. Results are shown for different choices of the renormalization and factorization scales. 
We use MSTW2008 parton distribution functions [2^]. Kinematic cuts on leptons, jets and missing energy are specified in 
Sec. IIVI The first column corresponds to purely vector Z' , the second one to purely axial Z 1 , and the third one is the reference 
point used in experimental analysis by the ATLAS and CMS collaborations [l7l. [22l|. 



LHC experiments to present their exclusion limits. 

In the context of these leptophobic top-color models, 
some of parameters are correlated. Indeed, the Z' cou- 
pling to SM particles is described by the Lagrangian 



C = ~oi cot I 



Oh Z 



(1) 



where g\ is the weak coupling constant and cot Or is 
the mixing parameter that encodes the overall devia- 
tion of the Z' couplings to quarks from the electroweak 
ones. The left- and right-handed currents are defined 
as j£ = i L j fl ,t L + b L ~f^b L - ulj^ul - d~Ll^d L , and 
Jr = h {tRlntR - u R j tl u R ) + f 2 (bRj^bR - dRjudR). In 
the top-color model, one would require /i > and/or 
/2 < in the right-handed current for top-color tilt- 

tag 123 . _ 

Exclusion limits are often presented assuming fx — 1 
and ji = which corresponds to the vector coupling of 
Z' to top quarks but, of course, other choices of these 
parameters lead to richer physics. For this reason, we 
consider four different coupling choices in what follows, 
(A, /a) = (1,1), (-1,-1), (1,0), (0,-1). The first and 
the second cases correspond to pure vector or pure axial 
Z' couplings; the remaining choices describe mixed cases. 

We can use the interaction Lagrangian Eq. ([1]) to com- 
pute the total decay width of the Z' boson. At leading 
order we obtain 



1 z 1 



a cot 2 Onrriz' ' 



8 cos 2 0w 
(! + /?)-(! 



(3 + f'l + 2/|) 



fi -6/0- 



(2) 



where /3 t = yl — 4m 2 /m|,. We therefore take fx, ji, 
r^r and the mass of the Z' resonance as free parameters 
of the theory and derive the overall coupling strength of 
Z' to quarks from them. We note that this leads to a 
strong correlation between the strength of couplings to 



quarks and the width of the resonance which may be 
less pronounced in other models. We will return to this 
issue later in the paper when discussing the importance of 
the interference between signal and background processes 
that occurs at NLO QCD. 

As we pointed out already, one of the goals of this pa- 
per is to present a realistic description of the hadronic 
Z 1 — > tt production process that combines QCD radia- 
tive corrections and top quark decays. To this end, we 
have implemented all the relevant top quark decay chan- 
nels. For simplicity, we will only discuss results for fully 
leptonic channels in what follows. We use the narrow 
width approximation for top quark pair production as 
described in Ref. [25j . The parametric accuracy of this 
approximation is 0(T t /mt) and its practical reliability 
is superb as long as the on-shell intermediate tt state is 
kinematically allowed [26j | . 

For a complete NLO QCD analysis of the pp —> Z' —> 
tt — > XX process in the narrow width approximation for 
t and t, we need the following ingredients: 



• NLO QCD corrections to pp 
ized t and t. 



Z' — > tt for polar- 



• NLO QCD corrections to the Z' decay width. 
Those corrections are potentially important to 
properly describe shapes of the tt invariant mass 
distribution when 0(a s ) corrections are taken into 
account. These corrections can be taken from the 
computation of 0(a s ) corrections to Z — > bb decay 
for massive b quarks reported in Ref. [27]. 

• NLO QCD corrections to top quark decays, includ- 
ing NLO QCD corrections to the top quark width 
and the W width. All these corrections were al- 
ready implemented in the calculation of pp — > tt 
reported in Ref. [25[ and we take them from there. 

• NLO QCD corrections to the QCD background pro- 
cess pp — ¥ ti. We use the implementation described 
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Narrow Z 1 , F^, ' 


= 15 GeV 


Broad Z 1 , Y 


LO) _ 

Z' 


= 150 GeV 




(1,1) 


(-1,-1) 


(1,0) 


(o,-i) 


(1,1) 


(-1,-1) 


(1,0) 


(o,-i) 


LO 


3.70 


3.53 


5.39 


2.00 


36.8 


34.3 


53.4 


19.6 


NLO production, qq 


4.93 


4.74 


7.17 


2.67 


49.4 


46.3 


71.6 


26.4 


NLO production, qg 


-0.15 


-0.14 


-0.22 


-0.08 


-1.4 


-1.2 


-2.1 


-0.7 


NLO production, gq 


-0.02 


-0.02 


-0.03 


-0.01 


-0.2 


-0.2 


-0.3 


-0.1 


NLO decay 


-0.60 


-0.57 


-0.87 


-0.31 


-6.1 


-5.6 


-8.8 


-3.2 


Interference with QCD 


~ 


-0.13 


0.01 


-0.06 


~ 


-1.3 


0.1 


-0.6 


Full NLO 


4.17 


3.89 


6.06 


2.20 


41.7 


38.1 


60.5 


21.8 


K— factor 


1.13 


1.10 


1.12 


1.10 


1.13 


1.11 


1.13 


1.11 



TABLE II: Cross sections (in fb) for pp — > Z' — 5* ti — > llvvjj + X at the 14 TeV LHC split into different contributions, for 
the reference scale = m Z ' — 1.5 TeV. We use MSTW2008 parton distribution functions. Kinematic cuts on leptons, jets and 
missing energy are specified in Sec. IIVI The leading order process, computed with NLO parameters (two-loop a s , NLO PDFs 
and widths) is included in the "NLO production, qq" category. 



in Ref. [25J. 

• Interference of the background qq — > ti and signal 
— !• Z — > ti processes that appears at NLO QCD 
for the first time. 



III. NLO QCD CORRECTIONS TO pp ->■ Z' 



a 



In this section, we describe some details pertinent to 
the computation of NLO QCD corrections to pp — s- Z' — » 
ti. Through 0(a s ), we can express the qq —> Z' — s- ti 
amplitude as 



2N C 



+ I SiqitSitlq — —SiqUjbitit ) ^2 



Nr. 



(3) 



where N c — 3 is the number of colors. The one-loop con- 
tribution to the amplitude Az' is given by the two color- 
stripped primitive amplitudes C\ and C2. C\ describes 
the amplitude with the same color flow as the tree-level 
one Afr ■ ; it is computed from vertex corrections to the 
production and decay stage of the process qq Z' — >• ti. 
It also receives contributions from the required renormal- 
ization constants. The function C2 receives contributions 
from box diagrams where a gluon and a Z' connect the 
initial qq and the final ti states. 

It is clear from Eq.Q that the amplitude C2 corre- 
sponds to the color-octet flow in the s-channel while the 
amplitude C\ and the tree-level amplitude A%) corre- 
spond to the color-singlet flow in the s-channel. For this 
reason, C2 does not contribute to the production cross 
section of Z' if we restrict ourselves to the underlying 
process qq — » Z' — » ti. On the other hand, we note that 
the color-octet amplitude can interfere with the back- 
ground process qq — > g* — > ti. Such interference is not 
expected to be large for the narrow resonance but it is not 
clear a priory if it is negligible for broader resonances. 



To be specific, consider the pure QCD scattering am- 
plitude for qq — > ti through one loop. Its color decompo- 
sition reads [25| 



Aqcd = 9 2 s (Si q tAti 9 ~ j^-S iq i 9 6i t i^j A 



(o) 

QCD 



9 S 



(4) 



where tree- and one- loop contributions are separated. It 
follows from Eqs. (|3I4[) that the tree-level QCD amplitude 

Aqcd can interfere with C2. We can also see from Eq. @ 
that, at one loop, .Aqcd contains the amplitude B2 whose 
color flow is equivalent to the color-singlet exchange in 
the s-channel. For this reason, B2 interferes with the 
tree-level amplitude for the Z' production A^) , but B2 
does not contribute to the pure QCD process qq —> ti 
due to the color structure of the qq — > g* — >• ti amplitude 
at leading order. Therefore if we consider the physical 
process qq — » ti that can be mediated by both gluon and 
Z' exchanges and if we allow for the interference effects 
between the two processes-that a priori are not obvi- 
ously negligible for broad Z' resonances-two amplitudes 
beyond what has been considered in the literature so far 
are required. We compute those amplitudes and include 
the effects of the interference between signal and back- 
ground processes in the discussion below. 

We note that we neglect contributions of the one-loop 
"anomaly-like" diagrams gg — » Z' — > ti that are medi- 
ated by quark loops. Such contributions are expected to 
be negligible for several reasons. First, the gluon flux at 
center-of-mass energies relevant for Z' production is ex- 
pected to be about one-tenth of the qq flux. Second, the 
anomaly-free nature of the effective Lagrangian Eq. ((T|) 
ensures that gg — > Z' vanishes identically if all quarks 
are massless. The quality of the massless approximation 
for top quarks in the gg — > Z' amplitude is controlled by 
the ratio of the top quark mass to the Z 1 mass which is 
quite small for values of the Z' mass that are of interest 
to us. We conclude that the double suppression caused 
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FIG. 1: Top-antitop invariant mass distribution in the reaction pp — » ti for m' z — 1.5 TeV at the 14 TeV LHC at leading and 
next-to-leading order in perturbative QCD. Left panel: narrow Z', T' z — 15 GeV and /i = l,/ 2 = 0. Right panel: broad Z' , 
T' z = 150 GeV and /i = — 1, f-2 = — 1. We use MSTW2008 parton distribution functions. Kinematic cuts on leptons, jets and 
missing energy are specified in Sec. IIVI The JGfactors are computed for the central scale choice fi — m z >- 



by the smallness of the gluon flux and the validity of 
the massless approximation for top quarks allows us to 
neglect the "anomaly" contribution. 

We also neglect the contribution coming from exchange 
of the standard model Z and photon in the qq annihila- 
tion channel, as they are expected to be tiny [28| in the 
invariant mass region m t j mz> ~ 1 TeV which is dis- 
cussed in this paper. We note that if required by some 
particular study, the inclusion of such interference con- 
tributions in our computation is straightforward since we 
work at the amplitude level and all coupling and masses 
that we use are free parameters. 

As a last remark, we would like to elaborate on the 
QCD corrections to the Z' width. In general, we write 
the Z' propagator in the unitary gauge as 

V t = p 2 ~T~- f~ [-Or + ^r) ■ © 

p k z — m z , + vmz'i- z 1 \ m z' / 

The longitudinal part of the Z' propagator is not im- 
portant since it always couples to the massless quark 
line on the "production side" of the process. When 
NLO QCD corrections to the production of Z' are con- 
sidered, it is essential to include corrections to T' z in 
the denominator of the Z' propagator of Eq. ([5]). To 
see this, consider production of a narrow Z', such that 
a(pp — y Z 1 — y ££) ~ a(pp -> Z') x Br(Z' ->• ti). To 
correctly compute the change in the branching ratio of 
Z' — > ti at NLO QCD, we need to account for correc- 
tions to the partial width Z' — > ti that we do by com- 
puting QCD corrections as described above and for the 
corrections to the total width Yz< that we have to com- 
pute separately and include in a Breit-Wigner propaga- 
tor. We read off the NLO QCD corrections to the Z' 
width from the calculation of QCD corrections to the de- 
cay rate Z — »• bb described in Ref. |27|, where they are 
given as functions of mb/mz- For numerical examples 
that we consider below mt/mz' <C 1, & n d in that limit 
corrections to the width become simple. Indeed, if we 
use the massless top quark limit, the NLO width is just 
(1 +a s /7r) times the LO width. Nevertheless, in our cal- 



culation, we retain full mt/mz> dependence of the width, 
following Ref. [27[ . 



IV. NUMERICAL RESULTS 

In this section we illustrate our computation by pre- 
senting some numerical results. We consider the LHC 
with yfs — 14 TeV center-of-mass energy. We consider 
semileptonic decays of both t and i. We use the anti-fct 
jet algorithm with R = 0.5 and impose cuts on the final 
state particles that are motivated by current Z' searches 
[l7L [22j |. In particular, we require p±,i e p < 20 GeV, 
|?7i ep | < 2.5, p±j Ct < 30 GeV, |jy jet | < 2.5. We choose 
mz' = 1.5 TeV, which is close to the current exclusion 
limits for our choices of the Z' couplings. We use the 
MSTW2008 parton distribution functions [H]. 

As pointed out in the Introduction, instead of spec- 
ifying the mixing angle 6h we find it more convenient 
to use the leading order width as an input parameter. 
We choose two reference values T z 9 = 15 GeV (narrow 
Z') and T z ? = 150 GeV (broad Z'). We also choose 
m Z ' as a reference value for the renormalization and 
factorization scales and show results for three choices 
(X = mz>/2,mz' and 2mz>- We note that, most likely, 
this procedure overestimates the theoretical uncertainty 
since it implies a scan over a very broad range of scales, 
between /x = 750 GeV and 3 TeV. 

We summarize our results for production cross sections 
in the dilepton channel for both narrow and broad Z' res- 
onances in Table HI The inclusion of NLO QCD correc- 
tions leads to extremely stable results under changes of 
the renormalization and factorization scales. The resid- 
ual scale dependence is of the order of 1%. The mild 
scale dependence is in part due to the fact that the LO 
result does not explicitly depend on a s so that the scale 
variation only comes from parton distribution functions. 
For fi = m Z /-we find that the NLO QCD corrections en- 
hance the production cross section by about 10%, so the 
corresponding .fG factor is 1.1. This value is somewhat 
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FIG. 2: Distribution of the effective mass defined as the sum of the transverse missing energy and H± for broad fi = f 2 = —1 
Z' at the 14 TeV LHC. We use MSTW2008 parton distribution functions. Kinematic cuts on leptons, jets and missing energy 
are specified in Sec. IIVI 



lower than K w 1.3, which was reported in Ref. [2J, and 
we attribute this difference to our use of MSTW2008 par- 
ton distribution functions, kinematic cuts on top quark 
decay products and inclusion of radiative corrections to 
top quark decays. We note that we can reproduce the 
results for the NLO QCD ii'-factors reported in Ref. @ 
if we use the setup for the calculation described in that 
reference. 

To understand the importance of various NLO QCD 
effects, we display in Table ITU all the different contribu- 
tions to the NLO cross section for n = mz'- The "pro- 
duction" contribution corresponds to NLO QCD correc- 
tions to various production processes pp — > Z' — > tt and 
leading order top quark decays. We include the leading 
order process qq — > Z' — > tt computed with the QCD- 
corrected values of T z> and T t in the production contri- 
bution. The "decay" contribution corresponds to leading 
order qq -> Z' -> tt and NLO QCD decays of t and t. The 
negative correction from the decay somewhat compen- 
sates the positive change in the leading order processes 
caused by using r^ LO there, but the cancellation is not 
complete because of the fiducial volume cuts. Finally, 
the last contribution displayed in Table HI] refers to the 
interference between QCD amplitudes and Z' amplitudes 
that arises at NLO QCD as we already discussed. 

An interesting feature of these results is that although 
the qq channel gives the largest contribution and it has 
an associated K- factor of about 1.3, similar to the result 
adopted in a recent analysis by the ATLAS and CMS 
collaborations fl7l . l22j ] , all other contributions to the K- 
factor are negative. 1 The net effect of these corrections is 
the reduction of the if-factor to a smaller value, K ~ 1.1. 

We also note that the relevance of the interference be- 



1 We mention that ATLAS and CMS use these if-factors for 7 
and 8 TeV LHC, whereas we have only shown results for 14 TeV 
LHC. We have, however, calculated the _fC-factor for 7 TeV LHC 
for the narrow Z 1 with a mass of 1.5 TeV for the coupling choice 
/i = l,/2 =0 and fiducial cuts as described in the text; we 
found K = 1.1 for fi = m z i. 



tween signal and background production processes de- 
pends on the details of the model. For the Z 1 model that 
we employ in this paper, the interference can reach a few 
percent, but it can be larger in general. The interference 
is absent for pure vector couplings of Z 1 due to Furry's 
theorem, and it is maximal for pure axial couplings. For 
the model specified by the Lagrangian Eq. ([T]) the mag- 
nitude of interference contributions can be as large as the 
production contribution in the qg channel. 

The interference effects become more relevant for 
broader Z 1 resonances. While this is an expected fea- 
ture of the result, it is perhaps not as apparent in Ta- 
ble [IT] as it should be, so we comment on how to prop- 
erly interpret what is seen there. Indeed, a glance at 
Table |n] suggests that the relative importance of the in- 
terference is similar for narrow and broad Z' . However, 
this is an artifact of our model choice where the growth of 
the Z' width is correlated with the growth of the Z' cou- 
plings to quarks. Since the production cross sections for 
pp — > Z' — > tt is quadratic in those couplings, while the 
interference is linear, the fact that both the cross section 
and the interference increase by an order of magnitude 
for the broad Z' compared to narrow Z' suggests that 
the relative importance of the interference increases by 
a factor of 3. Hence, we conclude that such interference 
effects can be important for broad resonances, especially 
if the width of the resonance is only weakly correlated 
with its coupling to quarks. We note that this interfer- 
ence contribution can play a role in such observables as 
the top quark forward-backward asymmetry at the Teva- 
tron. Although a flavor-conserving Z' resonance in the 
s-channel is already ruled out as a possible explanation 
for the observed asymmetry, it is interesting to note that 
the NLO QCD interference contribution to the asymme- 
try in our reference model can be as large as 20% of the 
leading order asymmetry generated by Z' exchange. The 
size of the effect suggests that it is worthwhile to study 
it for more realistic models that explain the Tevatron 
asymmetry. We leave this for future work. 

We now proceed to the discussion of selected differ- 
ential distributions. We begin with the distribution in 
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FIG. 3: Distributions of <j>+ variable, introduced in [23|], at leading and next-to-leading order in QCD at 14 TeV LHC. Both 
signal and background are included in the \m t i — m Z ' \ < 100 GeV window. Left panel: narrow (1,0) Z' model. Right panel: 
broad (-1,-1) Z' . We use MSTW2008 parton distribution functions. Kinematic cuts on leptons, jets and missing energy are 
specified in Sec. IIVI 



the tt invariant mass. While the ti invariant mass is dif- 
ficult to reconstruct in the dilepton channel, it is still 
instructive to look at the generic features of this dis- 
tribution. Such distributions for narrow and broad Z' 
are shown in Fig. [T] where we also include the contri- 
bution from the pp — > ti background process at lead- 
ing and next-to-leading order in perturbative QCD. We 
employ the reference parameter choices (/i,/a) = (1,0) 
and (fx, fz) = (—1, —1), respectively for the narrow and 
broad Z' cases reported in Fig.[TJ The m t j invariant mass 
distribution shows expected features-the Z' peak, which 
is very prominent at leading order, becomes broader and 
less prominent at NLO QCD since a significant fraction 
of events is shifted to lower m t i values because of the 
final state radiation. As a result, the if-factor changes 
significantly around the peak position, for both narrow 
and broad Z' . This rapid change in the K- factor may 
be invisible for the narrow Z', given that the resolution 
of m t i is close to 100 GeV but it may become important 
for broader resonances. Indeed, as follows from the right 
panel of Fig. [T] the if-factor changes by almost a fac- 
tor of 2, from K = 1.8 at m tl = 1300 GeV to K = 0.8 
at mti — 1600 GeV despite the fact that the average 
K- factor is relatively modest. Note that the if- factor 
at lower invariant mass values-which are more accessible 
experimentally-is actually quite significant and it almost 
doubles the expected cross section there. As a further 
illustration of this point, we note that the peak of the 
NLO QCD rati distribution for a broad Z' shown in the 
right panel of Fig. Q] is located in a different bin when 
compared to the leading order m t i distribution. 

Our calculation allows us to compute any kinematic 
distribution that is relevant for Z 1 production through 
NLO QCD and, in particular, those that are defined us- 
ing momenta of top quark decay products. For illus- 
tration purposes, we show two such distributions. In 
Fig. [2] the distribution of the effective mass for events 
with and without (broad) Z 1 is shown. All NLO QCD 
effects that are described above are included in the cal- 
culation. Understanding the effective mass distribution 



is of paramount importance for the experimental anal- 
ysis, where it is used to extract information about the 
Z' mass in channels where direct reconstruction of m t j is 
not possible. It follows from the the right panel of Fig. [2] 
that the if -factor is a rapidly changing function of m e ff. 
This is a consequence of the fact that the m c s distribu- 
tion receives contributions from signal and background 
processes which have markedly different dependence on 
m e g and very different if-factors. 

As another example, we show a kinematic distribution 
that can be constructed if lepton momenta are accurately 
measured. As was pointed out in Ref. (23|, such distri- 
butions are useful for exploring the Lorentz structure of 
the couplings of the new resonance to top quarks. Since 
our NLO QCD computation includes the spin correla- 
tions exactly, we can extend the discussion of the relevant 
kinematic distributions by including NLO QCD correc- 
tions. 

We consider the kinematic distribution of the az- 
imuthal angle 0+, which is constructed in the following 
way [23|. We work in the center-of-mass frame 2 of t and 
i. In that frame the directions of the incoming protons 
are not back to back; we denote them as ni and n%. We 
define a three-vector n\2 = (n\ — n2)/\n\ — r?2 1 and form 
a coordinate system by taking the z axis to be the direc- 
tion of the top momentum n t , the x axis to be n\2 and 
the y axis to be ri t x n 12 . We now write the momenta of 
the lepton and the antilepton in that reference frame as 
Pil= Ei i (l , sin $ cos (fii f, sin $ sin <p l j, cos 9 t j) and define 
the angle 0+ as 4> + = cos -1 (co s ipi cos ipj — sin ipi sin (pf) . 
As was pointed out in Ref. [23J, the 4>+ distribution can 
be used to discriminate between vector and axial cou- 
plings of Z' to a tt pair. 

We show the results of the computation of the </>+ dis- 



2 Because of the missing energy in dilepton events, it is difficult, 
but not impossible, to reconstruct the tt center-of-mass frame. 
See the discussion of this point in Ref. [2^1 . 
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tribution in Fig. [3] where two scenarios (i) narrow Z' with 
fi = 1 and f 2 — and (ii) broad Z' with f 1 = f 2 = —l 
are shown. The distributions contain t and i from both 
the QCD background process and from the Z' produc- 
tion. To reduce the contribution of the background, we 
require that the invariant mass of the ti system satisfies 
the constraint \m t i — mz>\ < 100 GeV. With this cut, 
the signal-to-background ratio is roughly 0.5 for the nar- 
row resonance and is close to a factor of 2 for the broad 
resonance case. It follows from Fig. [3] that the shape of 
this distribution is fairly stable against NLO QCD ra- 
diative corrections although in case (ii) (right panel) the 
QCD corrections make the distribution more flat. Since 
a small admixture of a vector coupling will have a similar 
effect on that distribution, it is essential to accommodate 
radiative corrections for the precise measurement of Z'tt 
couplings. 

V. CONCLUSION 

In this paper, we discussed NLO QCD radiative cor- 
rections to the production of a Z' boson in the reaction 
pp — > Z' — > ti. In contrast to previous studies, our com- 
putation includes interference effects of the Z' produc- 
tion process and the QCD background process that may 
become relevant for broad resonances. In addition, we 
include radiative corrections to top quark decays, keep- 
ing all spin correlations intact. We find that the NLO 
QCD corrections increase leading order cross sections by 
a modest amount that is somewhat smaller than what has 



been adopted in recent studies by the ATLAS and CMS 
collaborations [13, HH . We also find that radiative cor- 
rections to the ti invariant mass distribution and to the 
effective mass distribution are significant and strongly 
dependent on m t i and m g, respectively. For broader 
resonances, they increase the number of events with a 
smaller invariant mass of the ti system making heavier 
resonances more accessible experimentally. 

Our calculation is useful for understanding the Lorentz 
structure of the Z' interaction with the ti pair, with high 
accuracy, once the resonance is discovered. We have illus- 
trated this by computing the distribution in the relative 
angles of leptons and antileptons, defined as suggested 
in Ref. [23|. This variable was found [23| to be a good 
discriminator between vector and axial couplings of Z' to 
top quarks; our computation allows one to check if this 
conclusion remains valid when NLO QCD corrections are 
included. We found that, in the case of a narrow Z', the 
effect of QCD radiative corrections is to rescale this dis- 
tribution by an overall factor without changing its shape 
while the situation with broader Z' is somewhat more 
complicated. At any rate, our calculation makes it possi- 
ble to analyze any Z' scenario with NLO QCD accuracy. 
We look forward to the confrontation of our computation 
with LHC data. 
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